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Effect of chemical heat release in 
a temporally evolving mixing layer 

By F. J. Higuera 1 AND R. D. Moser 2 


Two-dimensional numerical simulations of a temporally evolving mixing layer with 
an exothermic infinitely fast diffusion flame between two unmixed reactants have 
been carried out in the limit of zero Mach number to study the effect of the heat 
release on the early stages of the evolution of the flow. Attention has been directed 
to relatively large values of the oxidizer-to-fuel mass stoichiometric ratio typical of 
hydrocarbon flames, and initial vorticity distributions thicker than the temperature 
and species distributions have been chosen to mimic the situation at the outlet of a 
jet. The results show that, during the stages of the evolution covered by the present 
simulations, enhancement of combustion occurs by local stretching of the flame 
without much augmentation of its area. The rate of product generation depends 
strongly on the initial conditions, which suggests the possibility of controlling the 
combustion by acting on the flow. Rollup and vortex amalgamation still occur in 
these reacting flows but are very much affected by the production of new vorticity 
by baxoclinic torques. These torques lead to counter rotating vortex pairs around 
the flame and, more importantly, in thin layers of light fluid that leave the vicinity 
of the flame when the Kelvin-Helmholtz instability begins to develop. Propelled by 
the vortex pairs, these layers wind around, split on reaching high pressure regions, 
and originate new vortex pairs in a process that ends up building large-scale vortices 
with a vorticity distribution more complex than for a constant density fluid. 


1. Introduction 

Combustion typically involves large amounts of energy release that very often 
result in strong couplings of the molecular transport and chemical reactions with 
the motion of the gas where such processes occur. A fast chemical reaction between 
two initially unmixed reactants in a mixing layer is a prototype of many practi- 
cal situations. In this case reaction occurs in a very thin flame controlled by the 
diffusion of the reactants from the carrying streams. Furthermore, the diffusion 
fluxes reaching the flame depend on the nature of the flow through the action of 
large-scale mixing and by the stretching and wrinkling that the flow imposes on the 
flame. Conversely, the flow may be very much affected by the gas expansion and 
baroclinic torques due to the chemical heat release. The coupled problem rapidly 
becomes very complicated and much work is required to analyze the elements of 
the coupling separately. 

1 E.T.S. Ingenieros Aeronauticos, Pza. Cardenal Cisneros 3, 28040 Madrid, Spain 

2 NASA Ames Research Center 
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Experimental investigations of gas-phase reacting shear flows under conditions 
of low heat release, with no apparent coupling between heat release and the fluid 
mechanics, have been conducted by Mungal k Dimotakis (1984), Mungal, Herman- 
son k Dimotakis (1985) and Mungal & Frieler (1988) using the hydrogen-fluorine 
reaction, and by Masutani k Bowman (1986) using the nitric oxide-ozone reac- 
tion. The amount of mixing and product formation, the effects of the Reynolds 
number and chemical reaction rate, and the structure of the reacting layer prior 
to the mixing transition are some of the issues addressed in these investigations. 
Models of the mixing and chemical reaction under these conditions were proposed 
by Marble k Broadwell (1977) and by Broad well k Breidenthal (1982) and Broad- 
well k Mungal (1991). On the numerical side, Riley, Metcalfe k Orszag (1986) 
carried out three-dimensional simulations of a temporally-evolving constant density 
flow with a single-step temperature-independent chemical reaction and showed that 
their results agree with similarity theory and compare well with the experimental 
data of Mungal k Dimotakis. Givi, Jou k Metcalfe (1986) addressed the problem 
of local extinction of a single-step irreversible Arrhenius reaction by means of two- 
dimensional temporally evolving simulations, and Ghoniem k Givi (1988) used the 
same kinetics in a two-dimensional spatially evolving simulation. A simulation of 
the spatially evolving flow with an infinitely fast reaction was performed by Delhaye 
et al (1994), who also proposed a model based on an equation for the mean flame 
surface density and a description of the strained flame elements. 

The effect of heat release on the flow was investigated by Wallace (1981), using 
the nitric oxide-ozone reaction with adiabatic temperatures up to 700 K and, more 
thoroughly, by Hermanson k Dimotakis (1989). They used the hydrogen-fluorine 
reaction with adiabatic temperatures ranging from 486 to 1240 K and showed that 
the growth rate of the layer, the turbulent shear stresses, and the overall entrainment 
are reduced as a consequence of the heat release. They also found that large-scale 
structures persisted in their experiment, but the mean structure space scaled with 
the layer width decreased with increasing heat release, which led these authors to 
suggest that the mechanisms of large-structure coalescence are inhibited by heat 
release. A review of further experimental work can be found in Dimotakis (1991). 

McMurtry et al (1986) studied the coupling between chemical heat release and 
fluid dynamics by means of direct numerical simulations of a two-dimensional, tem- 
porally evolving mixing layer with a single-step, irreversible exothermic reaction 
independent of the temperature. Using a zero Mach number formulation, they 
found that the rate of product generation, the thickness of the mixing layer, and 
the vorticity at the center of the vortex structures all decrease with increasing rates 
of heat release. They analyzed the influence of the baroclinic torque, pointing out 
that it changes sign across the flame because the direction of the pressure gradi- 
ent in the high temperature region is roughly radially outward from the centers 
of the vortices, while the density gradient changes sign across the reaction front. 
As the flow evolves, this results in the generation of vorticity that alternately op- 
poses and enhances the original vorticity, leading to several local extremas and to 
a vorticity distribution more diffuse than for a constant density fluid. McMurtry, 
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Riley & Metcalfe (1989) confirmed these results in a three-dimensional simulation 
and found that the Reynolds stresses, the turbulent kinetic energy, and the three- 
dimensionality of the flow also decrease with increasing heat release, in agreement 
with the experimental results of Hermanson & Dimotakis. Similar results were 
obtained by Grinstein &; Kailasanath (1992) in a two-dimensional simulation of a 
compressible subsonic spatially evolving mixing layer. 

Extensive simulations of compressible and/or reacting mixing layers have been 
carried out by Planche & Reynolds (1992) using an exothermic Arrhenius kinetics 
and values of the equivalence ratio such that the flame is near the center of the 
layer. They found that this flow behaves as two independent colayers separated by 
an approximately plane flame. Mixing occurs only between fuel and combustion 
products on one side of the layer and between oxidizer and products on the other 
side, but no roll-up or pairing was seen. Instead, the authors propose that the 
modification of the mean flow by the structures of each colayer causes the decay of 
these structures and the emergence of new ones of larger scale, which results in a 
slow growth of the layer. 

The purpose of this work is to investigate in further detail the effect of chemical 
heat release on the first stages of the development of a plane mixing layer at zero 
Mach number. As in much previous work, we restrict ourselves to two-dimensional 
numerical simulations, which should provide valuable information since it is abun- 
dantly clear that these early stages are dominated by the dynamics of large-scale 
two-dimensional structures. The chemical reaction takes place between two unmixed 
species and the reaction time is supposed to be much shorter than the shortest time 
scale of the flow, so that the molecular mixing of the reactants is confined to a con- 
tinuous infinitely thin diffusion flame. This assumption leaves out many important 
phenomena related to ignition and extinction but, in its broad range of applicabil- 
ity, gives results independent of any specific reaction law and renders the numerical 
resolution of the reaction zone unnecessary. A further limitation comes from the 
temporally evolving nature of the computations, which leaves out the effects of 
entrainment asymmetry (see, e.g., Dimotakis 1986 for a discussion of entrainment) 
and feedback from the downstream flow. While both effects are very important for a 
real reacting mixing layer, we expect that the temporally evolving flow still contains 
the essential ingredients of the internal dynamics of the mixing layer. As for the 
results, while a drastic reduction of the efficiency of the large-scale mixing seems to 
be unavoidable for non-small Mach numbers or when the flame sits near the center 
of the mixing layer, we find that this is not necessarily the case for the relatively 
high values of the effective stoichiometric ratio (or equivalence ratio) and initial 
vorticity distributions typical of many practical combustion systems. In this case, 
the mean location of the flame is displaced toward the oxidizer side and the presence 
of low density fluid modifies the evolution of the vorticity without suppressing the 
roll-up and pairing processes typical of constant density mixing layers. 

2. Formulation and numerical method 

Consider a plane two-dimensional temporally evolving unconfined mixing layer 
between a gas stream of velocity U and temperature T Q carrying a mass fraction y 000 
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of an oxidizer and another gas stream of velocity —U and temperature T/ carrying 
a mass fraction yj & of a fuel. The two species react in an exothermic infinitely 
fast chemical reaction so that a diffusion flame exists in the mixing layer where an 
amount of heat Q is released per unit mass of burned fuel. The velocity U is much 
smaller than the speeds of sound in either stream so that compressibility effects 
are negligible, but the density and temperature of the gas change due to the heat 
release (and to the temperature difference between the two streams if T/ ^ T 0 ). 

In what follows the velocities and lengths are scaled with U and 

= 217/max ( du/dy ) 0 , the vorticity thickness of the initial mean flow, the temper- 

atures and densities with the fuel-stream values Tj and />/, and the pressure with 

P f u 2 - 

The mass and momentum conservation equations are 

^ + V-(pv) = 0, (1) 

j^(pv) + V • (pvv) = -Vp + -^ {V • (£Vv) + N} , (2) 

where Re = p fUS M / pj is the Reynolds number, p = p/pf depends only on the 

temperature: p = T a , 



and the normal viscous stress Re~ 1 (p v + pj 3)V • v is included in the pressure. 

In the limit of zero Mach number, and neglecting the variations of the mean 
molecular mass of the gas mixture, the equation of state is 


pT = 1 . (3) 

The velocity and pressure are L-periodic functions of the streamwise coordinate 
x (the choice of the period is discussed later) and satisfy 

dp d<f> . , (A , 

u —* ±1 , pv — > oc±4> , -q * T a ±~^ f° r V ) (4) 


where <t> = j Jq f? dx dy is the outflow due to the gas expansion, which is split 

between the two sides of the mixing layer by means of the factors a+ = 1/(1 + T^ 2 ) 
and a_ = a_(. — 1. These values result from matching the mixing layer to outer 
wave regions where acoustic effects would be important. 

Assuming that the Lewis numbers of the two species are equal to unity, the energy 
and species conservation equations are 


L(T) = 


L(yj) = -u> , 


L(i to) = sv> , 


(5a, b, c) 
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where 

d 1 

L (*) = + v • (p v *) - ■ (p V *) . (6) 

\jf and y Q are the mass fractions of fuel and oxidizer, pjUw/6 u is the mass of fuel 
consumed per unit volume per unit time, s is the mass of oxidizer required to burn 
the unit mass of fuel (a = M 0 u/Mf for the overall reaction F + vO —* products, 
where M/ and M 0 are the molecular masses of the two species), and c p and Pr 
are the specific heat at constant pressure and the Prandtl number, both assumed 
constant. 

As can be seen by combining (5b)-(5c) and (5a)-(5b), the mixture fraction and 
the modified enthalpy, ( SY f -Y a + 1)/(1 + S) and (T - T a + qY f )/( 1 - T 0 + q) 
respectively, satisfy transport equations similar to (5) but without reaction terms 
in the right hand sides (Williams 1985). Here Yf = yj/y/oo and Y 0 = y 0 /y 000 , S = 
sy/oo/yooo is the effective stoichiometric ratio, and q = Qyjoo/cpTf. In addition, 
both variables satisfy the same boundary conditions and, assuming that the initial 
conditions are also the same, we find that 



SY f -Y 0 + l T-T 0 +qY f _„, jN 

1 + 5 “ l — To + q - Z ^y^ 

( 7 ) 

satisfies 

HZ) = 0 , 

(8) 

and 

Z = 0 for y — ► oo 
Z — 1 for y — ► —oo 

( 9 ) 


Z(x,y,t) = Z( x + L,y,t). 

In the Burke-Schumann limit of diffusion-controlled flames eqs. (5a-c) yield w = 0 to 
leading order (see Williams 1985 for details). The precise form of the consumption 
rate need not by specified in this limit other than requiring that w = 0 implies 


Y f Y 0 = 0 , (10) 

which means that the two reactants do not coexist, being separated by an infinitely 
thin flame where both concentrations vanish. Eqs. (7) and (10) suffice then to de- 
termine the temperature and species concentrations in terms of the mixture fraction 
Z. At the flame the mixture fraction and temperature take the values 

Z.t = — and T ad = T 0 + 

and 

Y s = *1*2=1 , Y 0 = 0 , T = 1 + (T ad - 1)±§£ 

17 = 0 , Y 0 = T = T 0 + (1 — T 0 + q)Z 




(ii) 


for Z > Z s % 
for Z < Z s x 


( 12 ) 
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Notice that since Z and its first derivatives are continuous across the flame, the 
gradients of oxidizer and fuel concentrations and the gradients of temperature on 
either side of the flame verify 


dYo = _ s ay> 

dn + dn 


and 


dT 


dT 

^ dY f 

dn 

+ 

dn 

_ +9 ^r 


(13a, b) 


where n is the normal to the flame. 

The problem (l)-(4), (8), (9), and (12) is to be solved subject to the initial 
conditions discussed in the following section. 

For the numerical treatment, a second order finite volume method on a staggered 
grid is used for eqs. (1), (2), and (8). The time integration, which is also second 
order accurate, is explicit for the convection terms, implicit for the transport terms, 
and uses a time splitting for the pressure correction. This amounts to the standard 
decomposition (pv)^ n+1 ^ — (pv)^ = ^i(pv) — A t^6p, and (1) yields 


A tV 2 6p = + V • [(pv) (n > + ^(pv)] . (14) 

The only difference with respect to a constant density fluid is the presence of 
(dp/dtY n+i ^ in the right hand side of (14). Since p = p(Z ), by (3) and (12), 
this term can be computed by evaluating (8) at the time level (n + 1). However 
(pvZ)(" +1 ) in the convection term of this equation gives the unknown contribution 


At 


V ■ (Z( n+1) fo) 

( jrfdz + z ) 


along with others already known (Z (n+1) is already known from the solution of 
(8), which can be computed before solving (14)). Thus, when this (dp/ 
is carried into (14) and into the boundary condition (4) for u, an elliptic problem 
results for 6p which differs from the usual one for a liquid. This problem is solved 
iteratively at each time step. 

Eq. (12) implies dT/dZ = (1 — T 0 + q) — (q/SZ s t)H(Z — Z s t), where if(C) 
is the Heaviside step function. In the numerical simulations, the discontinuity 
of the temperature derivative at the flame was smoothed out replacing H(() by 
Hp(Q = [tanh(/?C) + l]/2 and then, to the same approximation, T(Z) = T a d - 
(q/f3SZ 3 t) ln2 1 / 2 + (dT jdZ l )dZ\ which becomes exact when fi -> oo. This 
approximation amounts to assigning a finite thickness to the flame, measured by 
yS” 1 . It was checked that the results do not depend on ft provided it is sufficiently 
large for the transport effects to dominate in the scale of this artificial thickness. 
Typically ft = 20 in the computations presented below. 


3. Initial conditions and linear stability 

The initial conditions are the superposition of a base flow independent of x and 
small perturbations proportional to the eigenfunctions of the most unstable mode 
and its subharmonic, as given by a linear stability analysis of the base flow. 
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FIGURE 1. Temperature, velocity, and density-weighted vorticity for: (a) first type 
of base flow, and ( b ) second type of base flow. 

Two kinds of base flows are considered. The first is the result of letting an initial 
state consisting of a halfspace with Z = 0 and another with Z — 1, each moving 
with opposite velocities parallel to their common boundary evolve for a finite time 
t- This leads to a self-similar flow in which the mixture fraction and the streamwise 
velocity are of the form Z = Z(t)) and u = U(rj) with tj = f/* 1 / 2 and f = f y pdy, 
and satisfy 


[F(Z)Z V ] + —r)Z v = 0 , (15a) 

[F(Z)U V ] + = 0 , (156) 

Z = 0, U = 1 for Tj — ► oo , (15c) 

z = 1 , U = -1 for rj — ► -00 , (15cf) 


where F(Z ) = T a ~ l with T(Z ) given by (12), rounded as commented at the end 
of the previous section. The extra condition Z(0) = Z $ % is used to fix the origin; it 
amoun ts to measuring transverse distances from the flame. The transverse velocity, 
0(lfy/Ret) y due to the gas expansion, and the pressure are computed afterward 
from the equations of motion and the boundary conditions (4). The nondimensional 
vorticity thickness of this flow is = 2{tj Re) 1 ! 2 j and the condition 

= 1 determines the value of t in the previous expressions. The resulting profiles 
for q = 40, 5 = 9, and Pr = 1 are plotted in Fig. la. Notice that the flame (at 
y — 0) is displaced toward the oxidizer side of the mixing layer. This is because the 
oxidizer and fuel diffusion fluxes must satisfy (13a) and, for the present moderately 
large 5, this is possible only if the flame is in a region of small fuel concentration. 
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The previous solutions give vorticity and mixture fraction profiles of similar thick- 
ness when Pr — 0(1). This, however, is not a realistic representation of situations 
in which one or the two streams have been flowing parallel to a splitter plate or 
channel wall before forming the mixing layer. Then the vorticity at a short distance 
downstream the edge of the plate spreads over a boundary layer whose thickness 
depends on the length of the solid wall, whereas the diffusion of the species and 
the heat begins only around the trailing edge of the plate and extends therefore 
to a thinner sublayer. Assuming that the flow is laminar upstream, the thickness 
of this sublayer is 0(e* / 3 ) in the nondimensional variables of the previous section, 
where e = x/Re and x is the distance to the edge of the plate (Goldstein 1930). The 
structure of the laminar mixing layer for Re x Re */ 2 (if only one of the fluids 
move) or for Re » x » Re ~ 1 ^ 4 (if the two fluids move) consists of a self-similar flow 
in this sublayer and one or two thicker layers of cold vortical flow. This solution can 
be computed by using asymptotic matching. Here, with a view to mimic the flow 
in the case of a stream of fuel discharging into stagnant oxidizer without having to 
construct a uniformly valid solution out of the different asymptotic expansions, we 
consider a base flow with velocity u = 1 — 2 df /d£, where /(() and Z(£) satisfy 


with 


e[F(Z)Z^ + -PrfZ' = 0 

+ 

z = h = o 

Z = h = 1 


? G+ l e( l- G ) 

for £ — > oo , 
for i — ♦ — oo , 


fh ~ 3 Gf , 2 = 0 



1 1 + tanh 


7 




(16a) 

(166) 

(16c) 

(16<i) 


1 e 1 / 3 and 7 constants, and the condition Z( 0) = Z s t to fix the origin. 

Eqs. (16a, b) reduce to the ones describing Goldstein’s self-similar flow in the dis- 
tinguished limit [/ = e 2 / 3 /„( = e 1/3 &, (/.,(«) = ^(1)] for e -► 0, in which G -♦ 1, 
and to Blasius’ equation for an isothermal flow (Z — 0) in the distinguished limit 
(/,£) = 0(1) for e -+ 0, in which G -* 0. Thus, in an approximate sense, the solu- 
tion of (16) for small values of e describes the whole flow. The artificial parameters 
— ( c and 7 control the position and abruptness of the transition between the two 
regions. Fig. lb shows the solution of (16) for q = 40, S = 9, Pr = 1 , c = 2.2x10 , 

( c = 0.33, and 7 = 1, normalized to have 6 U = 1. 

A linear stability analysis of both types of flows is carried out using the parallel 
flow approximation and neglecting transport effects for the perturbations. Squire’s 
theorem holds for these zero Mach number flows and, therefore, it suffices to consider 
planar perturbations of the form tt = u s (j/) + Au(x,y), etc., where the subscript 
s denotes here the base flow and A u(x,y) = u(y)exp [ia(r — ct)] <C 1, with a real 
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FIGURE 2. Growth rate as a function of the wave number for (a) the first type of 

base flow and ( b ) the second type of base flow. In ( a ), the cases depicted are , 

9 = 40, S = 9, cr = 0; , q = 40, S = 9, a = 1/2; , q = 20, 5 = 4, a = 0; 

and , q = 0. In (6) the cases depicted are , q = 40, 5 = 9, a = 0; and 

: q — 40, S = 9, cr = 1/2, and the phase speed is on the right-hand scale. 


and c complex. After standard manipulations of the linearized equations, we end 
up with the eigenvalue problem 


(«J - c) (p,v')' - p s a 2 v 
v 0 for y —¥ ±oo , 


= 0 , 


(17) 


where v is the transverse velocity perturbation. The resulting growth rate (ac;) as 
a function of the wave number is plotted in Fig. 2a for several base flows of the first 
type discussed above and in Fig. 2b for base flows of the second type. 

It was pointed out by Shin &: Ferziger (1991) and Planche & Reynolds (1992) 
that the density- weighted vorticity —p 3 u f s plays to a large extent the same role as 
the vorticity for a constant density flow. This quantity has two distinct peaks for 
sufficiently exothermic reactions with values of S not fax from one, leading to two 
different unstable modes for a range of a bounded away from zero. Each mode is 
associated with one of the peaks and has a phase velocity closer to the velocity of 
the stream nearer to that peak. A remnant of the two peaks can be seen in Fig. la, 
and two unstable modes still exist for q = 20, S = 4 (Fig. 2a). However, the lower 
peak dominates for higher values of S, leaving only one unstable mode for each 
wave number. The corresponding eigenfunctions have their extremes around the 
flame position and around the peak of the weighted vorticity. 

Another important feature is that the phase velocity of the perturbations depends 
on the wave number; it is negative for the cases of Fig. 2a and, contrarily to what 
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Case 

a (/i oc T°) 

A, 

A2 

¥>1 

Re 

Base Flow 

0 ) 

0.0 

0.05 

0.05 

0 

1000 

1 

(*) 

0.0 

0.15 

0.05 

0 

1000 

1 

(c) 

0.5 

0.05 

0.05 

0 

1000 

1 

(d) 

0.5 

0.15 

0.05 

0 

1000 

1 

(e) 

0.5 

0.15 

0.05 

7r/2 

1000 

1 

(/) 

0.5 

0.15 

0.05 

0 

2000 

1 

(0 

0.5 

0.05 

0.05 

0 

1000 

2 

(D) 

0.5 

0.15 

0.05 

0 

1000 

2 

Table 1. Conditions for the various cases run. 

A\ and A2 are 

the subharmonic and 

fundamental amplitudes (see 
Base flow types (1 and 2) are 

text) and <pi is the relative phase of the subharmonic. 
; described in §3. 


can be seen in Fig. 2b, its absolute value increases with a. In general, the maximum 
growth rate increases by decreasing q or increasing S. 

4. Numerical results and discussion 

In all the cases discussed in this section the temperature of the two streams are 
equal (T 0 = 1) and, unless otherwise specified, q — 40 and 5 = 9, leading to an 
adiabatic temperature of 5 times the free-stream temperature. This relatively high 
value of 5 is typical of many real situations; e.g., 5 ~ 15 for hydrocarbons burning 
in air. The Reynolds number is 1000 in all the computations except case (/) below, 
and Pr = 1. The length L of the computational domain is twice the wavelength 
of the most amplified small perturbation for the corresponding base flow. The 
initial perturbation is given by the eigenfunction of this mode normalized so that 
max |u| = 1 and multiplied by the amplitude A 2 = 0.05 plus that of its subharmonic 
taken with various amplitudes (A\) and phases (y?i). The specific conditions of each 
case are listed in the Table. The cross-stream boundary conditions are applied at 
y = ±8 in computations with the first type of initial conditions and at y = —9 and 
y = 7 in computations with the second type of initial conditions. 

4.1. Global results 

Global results for a number of simulations carried out with the first type of initial 
conditions axe displayed in Fig. 3. Fig. 3a shows the time evolution of the mean 
momentum thickness (6 m = /^^(l — 'pu 2 )dy, the overbar meaning streamwise 
average) seeded with its initial value. As can be seen, the growth rate of the mixing 
layer is similar for all the cases displayed, and smaller than for a constant density 
fluid. Fig. 3b shows the overall rate of product generation (W), 
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where 

v _ J (1 - Z)/(l - Z.t) if Z > Z,t] 
p \ ZfZ't if Z < Z tt ; 

is proportional to the product mass fraction. Figs. 3c and 3d show the first two 
Fourier modes of the transversely averaged kinetic energy 

Ek = \J j(? 1/2 “) (*> V) + (/> 1/2 v) (*,y))dy , 

where the hat denotes streamwise Fourier transform. The modes in Figs. 3c-d 
correspond to the wave numbers of the most amplified perturbation (Jb = 2) and its 
subharmonic (k = 1). 

The first peaks of the curves of Fig. 3c represent the roll-up of the vortices, which 
is accompanied by a first maximum of the product generation, and the first peaks 
of the curves of Fig. 3d represent the amalgamation of pairs of vortices. Cases (a) 
and (6) correspond to a constant viscosity fluid (<r = 0) with the amplitude of the 
subharmonic A\ = 0.05 and 0.15, respectively. The amalgamation occurs earlier 
and the reaction rate has a tall peak shortly after the maximum of E\ when the 
amplitude of the subharmonic is larger. Cases ( c ) and ( d ) are the analogs of (a) 
and ( b ) with a variable viscosity fluid ( a = 0.5), whereas case (e) differs from ( d ) 
only by a 7r/2 shift in the phase of the subharmonic. The effect of the initial phase 
shift is less obvious here than for a constant density fluid because the phase speed 
of the perturbations depends on the wave number and, therefore, the phase shift 
changes with time. However, a strong influence remains, as the change in phase 
results in a delay in the amalgamation and suppresses the tall peak of the rate of 
product generation. Case (/) has the same initial conditions as (<f) but Re = 2000. 
The differences between the two reflect the influence of the viscosity, which may be 
strong at Re = 1000, especially when the viscosity increases with temperature. 

In none of the cases presented does the flame sheet enter the cores of the vortices, 
and its surface never increases very much. This is because the relation (13a) implies 
that the flame must always be in regions of small fuel concentration since 5 is 
moderately large, and such regions are not easily ingested by the vortices when 
cushioned by a layer of low density fluid (see next subsection). Similarly, owing to 
(13b), the temperature gradient on the oxidizer side must be much larger than on 
the fuel side when q is large, a fact that is confirmed in Figs. 5 and 6 below. Under 
these conditions, the peaks along the flame in the product generation rate axe due 
to the enhancement of the diffusion fluxes in regions of high strain on the flame 
surface. The strong dependence of this strain (and therefore the rate of product 
generation) on the amplitude and phase of the initial perturbations suggests that 
it is possible to control the combustion by manipulating the flow. 

Further computations were carried out decreasing q (hence the flame tempera- 
ture) at constant 5, and decreasing q and S simultaneously to bring the flame closer 
to the center of the mixing layer while keeping the flame temperature constant. The 
results of the first series of computations, similar to (a) but with q equal to 20 and 
30, show features approaching those of an incompressible flow as q decreases: the 
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FIGURE 3. Evolution of (a) the scaled momentum thickness, (6) the overall rate of 
product generation, (c) the kinetic energy in the second Fourier mode, and (c?) the 
kinetic energy in the first Fourier mode, for cases: (a); (fc); (c); 

— - (< 0 ; («); (/); «■ 


vortices get rounder than in Fig. 5 below, the layer grows faster, and the fuel lean 
region containing the flame is more easily wrinkled and stretched by the flow, which 
results in higher rates of product generation. The second series of computations, 
(q,S) = (20, 4), (12, 2), and (8, 1), show that the amalgamation is retarded or sup- 
pressed and the rate of product generation decays with oscillations after peaking at 
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FIGURE 4. Evolution of (a) the scaled momentum thickness, ( b ) the overall rate of 
product generation, (c) the kinetic energy in the second Fourier mode, and ( d ) the 

kinetic energy in the first Fourier mode, for cases: (C); (D); (d); 

(d'); and a constant density flow. 


roll-up when 5 begins to decrease (case (g) in Fig. 3b is the first case of this series; 
for still smaller values of S the roll-up also disappears, apparently giving way to the 
two-colayers regime discussed by Planche &; Reynolds (1992)). 

Fig. 4 shows the same quantities as Fig. 3 for two runs carried out with initial 
conditions of the second type: cases (C) and (D) are analogous to (c) and (d); i.e., 
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a — 1/2 and the amplitude of the subharmonic equal to 0.05 and 0.15, respectively. 
For comparison, the results of case ( d ) have been repeated in this figure (dashed 
lines), and the results for a constant density fluid (q = 0, 5 = 9) have also been 
included (dotted lines). Finally the dash-and-dot lines correspond to case ( d ') dis- 
cussed below. As can be seen, the growth of the mixing layer is faster than for 
the first type of initial conditions, and similar to that of a constant density fluid 
for times up to about 20. The amplitude of the second mode of the kinetic energy 
(E 2 ) is also higher, of the order of the amplitude of the first mode (£ 1 ), reflecting a 
more efficient roll-up and the presence of round vortices at least during part of the 
evolution. The two peaks of E 2 reflect smother characteristic feature of this flow 
that will be commented on below. 

4.2. Mechanics of the flow 

Fig. 5 shows the evolution of the vorticity and the temperature for case (/), which 
is typical of the first type of initial conditions. In accordance with the results of the 
linear analysis, the vortices move with a velocity closer to that of the fuel stream 
(the lower stream in Fig. 5) than to that of the oxidizer stream (the structures move 
to the left). Furthermore, the interaction of the vortical region of the layer with 
the fuel side (lower) is typical of the constant density flow, with large excursions 
of free-stream fluid into the layer, whereas on the oxidizer side (where the flame 
is), these large-scale incursions are not so pronounced. The hydrodynamic mixing 
layer grows by entraining free-stream fluid as the vortices roll-up and amalgamate. 
On the oxidizer side, however, the heat release of the flame and the associated 
volumetric expansion of the surrounding fluid acts as a source of fluid volume to be 
entrained. Thus this low-density fluid below the flame, which is easier to entrain 
than the free-stream cold fluid, acts as a buffer between the mixing layer and the 
flame, inhibiting the entrainment of the flame and free-stream fluid. 

The light fluid enters each vortex in the form of a jet to its upper right, turns 
around the vortex to the left on the lower side, and splits into a recirculating flow 
around the vortex core and a blob of light fluid that leaves the vortex to the lower 
left (t = 9 to 15). The heavy fluid entering a vortex at its lower left turns over the 
top to the right and tends to escape as a tongue by its upper right (t — 15). The 
light fluid appears to be more easily turned and engulfed by the vortices than the 
heavy fluid, as might be expected given the difference in density. The right-ward 
moving tongue at the top throttles the ingestion of new light fluid from the upper 
right ( t = 15), and similarly, the left- ward moving blob of light fluid temporarily 
suspends the ingestion of new heavy fluid. 

The general effect of these combined motions of light and heavy fluid is to make 
the temperature field structures more elongated than for a constant density fluid. 
The numerical results show a variety of situations depending on the order in which 
the supplies of light and heavy fluid to a given vortex are interrupted and restored, 
which in turn depends on the initial conditions and on the Reynolds number. The 
supply of light fluid can be interrupted more than once before two neighboring 
vortices amalgamate. Thus, in case (6) (not displayed), a tongue of heavy fluid 
separates completely from a vortex, restoring the ingestion of light fluid, but it 
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remains cold on leaping over the next vortex to the right (which is undergoing a 
similar evolution), so that, along with some newly ingested fluid from the lower 
stream, it throttles the supply of light fluid to this vortex for a second time. For 
their part, the blobs of light fluid in the lower part of the layer cool down while 
moving toward the left and undergo a further splitting in this case, being partially 
ingested by the next vortex to the left of that which shed the blob. 

The right-ward motion of the tongues of heavy fluid over the vortices also deter- 
mines the end of the roll-up process, and when this happens sections of the original 
vorticity layer are left unrolled between the tongues and the flame. It seems likely 
that through this mechanism the phase of the subharmonic in the initial condi- 
tion may induce differences in the amount of vorticity rolled in each vortex and, 
thereby, in the subsequent evolution of the mixing layer. The unrolled vorticity 
patches correspond to regions of strong shear that move toward the right relative 
to the vortices, leap from vortex to vortex, and are selectively reinforced by the 
baxoclinic torque. At sufficiently high Reynolds numbers some of these patches 
eventually roll up into secondary vortices located one or more periods to the right 
of the vortex to which they were initially tied (upper left of Fig. 5b at t = 18 to 
21 ). 

The evolution leading to the premature termination of the roll-up and the dis- 
ruption of the rolled and unrolled vorticity can also be described in terms of the 
vorticity; see comments on Fig. 6 below. Here we note that in general the dis- 
tributions of both V x v and V * v are required to reconstruct the velocity field. 
The vorticity is initially present in the flow, and its dynamics are essentially in- 
viscid through convection and the baroclinic torque. In contrast, in the absence 
of compressibility effects, the changes in the density of a fluid particle are due to 
temperature changes, which are brought about by heat conduction. Thus at suffi- 
ciently high Reynolds numbers, heat conduction, and hence V • v = , should 

be confined to very thin layers (which can only have a small effect on the large-scale 
dynamics of the flow) and, in a fully developed flow, to the smallest scales (which 
cannot influence the large-scale dynamics by this means because the net change of 
volume is negligible). This points toward a decreasing importance of V • v as the 
Reynolds number increases. Though Re = 2000 is probably not sufficiently high 
to claim independence of V • v, the evolution depicted in Fig. 5 is probably more 
affected by the motion of the counterrotating vortex pairs that appear in the flow 
than by the dilatation or contraction of the fluid particles. Such pairs arise due to 
the baroclinic torque (— VT x Vp in the vorticity equation), which, owing to the low 
pressure at the vortex centers, is positive in the lower halves of the light fluid jets 
entering the vortices find in most of the heated region above the flame, and negative 
in the upper halves of the jets and below the flame. This leads to production of 
positive vorticity in the former regions, which forms left- ward moving pairs with 
the rolled vorticity of the layer and, to a lesser extent, right-ward moving pairs with 
the unrolled vorticity. 

It seems that in some cases, though not in the one displayed in Fig. 5, the unrolled 
layers of vorticity act as precursors to the amalgamation by establishing bridges 



Figure 5(a). Evolution of the temperature in case (/). The contour increment 
is 0.5 and the free-stream temperature is 1. The flame location is indicated by the 
heavy dashed line. 


between adjacent vortices with the right end of a bridge rolling around the vortex 
that will come to the bottom of the pair while the left end is still connected to 
the vortex that will come to the top. Any nearby secondary vortex that might 
have formed as well as the remnants of warm fluid blobs shed by other vortices are 
typically ingested in the amalgamation process. 

The result of the amalgamation is a vortex much less organized than for a constant 
density fluid. It contains several vorticity peaks and layers of positive vorticity 
left over by the secondary vortices and vortex pairs that existed at earlier times 
(t = 27 to 30). However, the large-scale behavior of the new vortex does not differ 
qualitatively from that of the previous ones: as far as it can be followed in the 
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FIGURE 5(6). Evolution of the vorticity in case (/). The contour increment is 
0.6 and positive contours are dotted. The flame location is indicated by the heavy 
dashed line. 


present simulations, the vortex grows by ingesting light and heavy fluid until it 
ends up filling the computational domain. 

To try to ascertain the importance of the baroclinic torque, the momentum 
Eq. (2) was artificially modified, replacing the pressure gradient term by — /?Vp. 
This amounts to having a constant density fluid insofar as the pressure forces are 
concerned. No baroclinic term appears then in the vorticity equation (obtained by 
taking the curl of the modified momentum equation), but the term -wV • v due 
to the gas expansion is left unaltered, as well as the effect of V • v in the other 
conservation equations. Case (d) was rerun with this modification (it is denoted by 
(d f ) in what follows) and some results are shown in Fig. 4 (dash-and-dot lines; the 




FIGURE 6(t). Evolution of the vorticity in case ( D ). The contour increment is 
0.8 and positive contours are dotted. The flame location is indicated by the heavy 
dashed line. 
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dashed lines correspond to case ( d )). As can be seen, the growth rate of the layer 
and the rate of product formation increase very much due to the modification, and 
inspection of the vorticity fields reveals that the vortices and the general appearance 
of the flow are closer to a constant density fluid than to case (d). 

Fig. 6 shows the evolution of the vorticity and the temperature for case (D), 
the analog of (d) with an initial condition of the second type. Contrary to Fig. 5, 
the bulk of the initial vorticity is now in the low temperature fuel. The Kelvin- 
Helmholtz instability begins to develop as for a constant density fluid, leading to 
relatively round structures (for t = 9 and before) that are responsible for the first 
peak of E<i in Fig. 4c. Very soon, however, the deflection of the light fluid around 
the flame leads to tongues where new vorticity is generated by the baroclinic torques 
in the form of counterrotating pairs. It is worth noting that this occurs in regions 
that would be relatively free of vorticity in a constant density fluid. The vortex 
pairs are more prominent than in Fig. 5 and clearly dominate the flow. They travel 
toward the left deep in the fuel stream (t = 11.25) before the dominant negative 
vorticity brings them back toward the shear layer, and in so doing, they make the 
elongated structures responsible for the dip of E 2 and engulf large amounts of cold 
fuel that are subsequently incorporated into the mixing layer. On returning to the 
shear layer, each of the regions of light fluid has the appearance of a mushroom with 
a very thin stalk. The high pressure above the cap of the mushroom splits it into a 
recirculating tongue and a shed blob, each with its associated vortex pair. The high 
pressure also squeezes substantial parts of the high temperature strip around the 
flame, increasing the diffusion fluxes and the burning rate while the parts of the strip 
between two pressure maxima get compressed and sucked by the depressions at the 
centers of the forming vortices. This leads to secondary tongues of light fluid (that 
begin forming at about t = 13.5 in Fig. 6a) where the baroclinic torques generate 
new counterrotating vortex pairs. Between t = 15.75 and t = 18, each secondary 
tongue is pinched off by a high pressure region and the recirculating part of one 
of the primary tongues, which splits again in this process into a layer leaving the 
vortex by its upper right and another that winds around along with the segmented 
secondary tongue. The process of formation of a secondary tongue and splitting of 
the recirculating flow occurs once again during this computation, and since there 
are several vortex pairs winding around, the vortex cores in Fig. 6b have for some 
time a multilayered structure until diffusion cancels vorticity of opposite signs (and 
temperature differences) leaving cores of predominantly negative vorticity (t = 27). 

Finally these cores undergo a pairing, during which one or two more counterro- 
tating vortex pairs are formed and ingested from the upper right of the forming 
vortex. At later times (not shown) these new pairs will become the front of a new 
tongue that heads toward the lower core moving toward the lower left while light 
fluid comes in behind and rolls around the surviving upper vortex. The result is an 
elongated structure with an associated decrease of E\ in Fig. 4d. 

Parts of the high temperature strip around the flame are alternatively strained 
and compressed by the evolving flow, and patches of fluid separate from this strip 
to form the tongues, but as was mentioned before, the flame itself remains outside 
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the vortex cores at least in the part of the evolution covered by these computations. 
This fact makes the flow in the vortices below the flame similar to a mixing layer 
between two fluids of different densities. To assess the extent of the similarity, 
further computations were carried out for an inert flow: q = 0 in the formulation 
of sec. 2 but T 0 = 5, so that the maximum temperature of the fluid equals the 
flame temperature of the previous reacting flow simulation. Except for the absence 
of baroclinic vorticity production around the flame, the results of this simulation 
(not displayed) have the same general appearance as those of Fig. 6, though minor 
differences do occur. 
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